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I. INTRODUCTION 

In the previous part of this article [1 (Part I), we 
have analyzed the correspondence between classical and 
quantum signatures of chaos in the Geometric Collective 
Model (GCM) [2 of nuclear vibrations. Rotations were 
ruled out by the constraint of zero angular momentum. 
The classical version of this model was previously shown 
[3 to exhibit a very complex dependence of regular and 
chaotic measures on control parameters and energy. 

In Part I we have compared the classical measure 
/reg, a regular fraction of the phase space volume, with 
the adjunct {l—oo) of the Brody parameter. The analysis 
of spectra was performed in a wide energy domain and 
for several values of the control parameter B of the GCM 
Hamiltonian. Spectra obtained via different quantization 
schemes of the classical model were considered, which 
led to the use of three different sets of quantum levels, 
denoted as 5D, 2D even, and 2D odd. Whereas the 5D 
spectrum corresponds to the standard five-dimensional 
GCM restricted to the nonrotating case, the 2D even and 
2D odd spectra were derived from the quantization in the 
two-dimensional space of polar deformation coordinates 
P and 7, with the respective condition on the parity of 
wave functions under the reflection of angle 7. 

In all cases, the validity of the Bohigas conjecture [4] 
has been fully confirmed. We stressed two important as- 
pects of our calculation: First, Bohigas' conjecture has 
been verified independently of the quantization scheme. 
Second, because of the strong dependence of the GCM 
chaotic measures on energy, the competing types of level 
statistics have been analyzed in the local regime, i.e. for 
separate portions of the spectrum. Whether the statis- 
tics in a given portion is more of the Poisson or Wigner 
type depends on the character (regular or chaotic, re- 
spectively) of the classical dynamics in the correspond- 
ing energy interval. The nonmonotonous dependence of 
chaotic measures on energy is in contrast to a majority 
of systems used as case examples of chaos, in particular 
to all kinds of billiards (or cavities) for which the chaotic 
features are energy independent. 

In this part of the contribution, we continue the work 
initiated in Part I by considering more sophisticated tech- 



niques to describe chaos in quantum spectra. Given that 
the Brody parameter captures only the short-range spec- 
tral correlations, a natural way to extend the previous 
results would be to consider also some measures of the 
long-range correlations, like e.g. the A3 statistics or the 
number variance [F. This way we did not follow. The 
reason is connected with the above-mentioned nontrivial 
variation of chaotic measures with energy, which would 
unavoidably increase statistical ambiguity of such analy- 
ses. 

Instead, we employed the method invented in 1984 by 
Peres [6 . While Bohigas' conjecture, which was pub- 
lished in the same year [4 , has become a widely rec- 
ognized paradigm of quantum chaos, Peres' method has 
been more or less forgotten. It is mentioned in the text- 
book [7 , where some applications in integrable and non- 
integrable spin systems are discussed . An application 
in a billiard system was presented later [9] . Today, how- 
ever, Peres' name is more commonly cited in connection 
with his alternative definition of chaos in quantum sys- 
tems [10^, submitted and published with a difference of 
just few days, which became a cornerstone for presently 
a quickly expanding branch of the quantum information 
theory [11 . 

Nevertheless, the idea of Ref. [6 turns out to be 
very fruitful, as well. We will show below that the 
method based on this idea represents a sensitive probe 
into the competition between regular and chaotic fea- 
tures in quantum spectra. It can be applied even beyond 
the theory of quantum chaos, as a synoptical indicator 
of the changing structures across the spectrum contain- 
ing possibly a very large number of states. In the two- 
dimensional case, the method is graphical and may be 
compared to the classical method of Poincare sections. 
The spectrum of stationary states of a given quantum 
system with two degrees of freedom is drawn as a lat- 
tice in the plane E x (P), where E is energy and (P) 
stands for an arbitrary observable average. This allows 
one to recognize ordered and disordered patterns and vi- 
sually allocate regular and chaotic domains within the 
same energy interval. The freedom in choosing observ- 
able P makes it possible to focus on various properties of 
individual states and to closely follow the way how chaos 
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sets in and proliferates in the system. 

The plan of the present paper is the following: Peres' 
method is introduced in Sec. [TTl and elaborated for the 



geometric model in Sec. Ill Section IV shows Peres lat- 
tices calculated with two choices of operator P for differ- 
ent values of control parameters and different quantiza- 
tions. Discussion of these results and their comparison 
with results of Part I are presented. The summary and 
conclusions come in Sec. [Vl 



II. PERES METHOD 

Let us consider an integrable system with two degrees 
of freedom (a 2D system). Apart from the Hamiltonian 
i^o, there must exist another integral of motion, denote it 
/, which by definition satisfies the commutation relation 
[Hq^I] =0. If we plot the eigenvalues of observable / 
against energies Ei for individual levels (enumerated by 
integer z = 0,l,2,...), the resulting image is formed by a 
lattice of regularly distributed points. This is a straight- 
forward consequence of the semiclassical quantization by 
Einstein, Brillouin, and Keller (EBK) [12 . An example 
of such a regular lattice is shown in panel (a) of Fig. [l] 
with Hq and / described in Sec. Ill In the present section 
we focus on the left column of the figure until specified 
otherwise. 

The kind of lattice described above is commonly used 
to analyze some analytic aspects of quantum integrable 
systems, see e.g. Ref. [T^. In the present context, it is 
a natural starting point for the explanation of the Peres 
method. To continue, we introduce a perturbation H' to 
the integrable Hamiltonian Hq, which yields a noninte- 
grable Hamiltonian 



(1) 



with A standing for a real number measuring the strength 
of the perturbation. Obviously, / does not commute with 
H since [H^I] 7^ in general, hence / is not any more 
an integral of motions. Consequently, the Hamiltonian 
eigenstates {i/ji) are characterized by energies Ei, but not 
by fixed values li of observable /. If we want to continue 
with the 2D visualization of the spectrum as introduced 
above for A = 0, the question raises what to draw on 
the vertical axis instead of J^? Quite naturally, Peres has 
chosen the expectation values = of observ- 

able / in individual eigenstates. This choice smoothly 
connects the perturbed (A 7^ 0) and the unperturbed 
(A = 0) cases since {I)i li for A ^ 0. 

Such lattices are shown in panels (b)-(e) of Fig. [l]( left), 
where the general parameter A from Eq. ([T]) is replaced 
by the model-specific parameter see Sec. |III[ As can 
be seen in panels (b)-(d), adding a small perturbation 
to the integrable Hamiltonian does not instantaneously 
break down the entire regular lattice. Instead, some lo- 
calized seeds of distortion are created, while the rest of 
the lattice remains ordered in the same fashion as in the 



integrable case. This scenario is in accordance with Per- 
cival's conjecture [M] assuming that the sets of regular 
and chaotic eigenstates are statistically independent in 
the semiclassical limit ^ 0, i.e. they do not interact 
with each other. Therefore, the persisting regular parts 
of the lattice can be associated with surviving remnants 
of classical tori, while the disordered parts correspond to 
proliferating chaotic orbits. 

As the perturbation strength A grows, the remnants 
of tori are gradually disappearing and disorder tends to 
increasingly plague the lattice. This is demonstrated by 
an almost totally disordered lattice in panel (e) of Fig. [ij 
where only a few low lying states keeps the regular pat- 
tern. 

The above-outlined visual method implies a great 
heuristic gain. It allows one to judge which parts of a 
mixed spectrum (or, in optimal cases, which individual 
states) are regular and which are chaotic. Let us stress 
that this is opposite to traditional methods of quantum 
chaos based on the spectral statistics since in that case 
regular and chaotic (or mixed) parts of the spectrum can 
only be specified by energy. In the present approach, 
these parts can coexist within the same energy interval, 
the additional information needed for their separation 
being obtained from the behavior of the averages (/)^. 

In practice, there certainly exist severe limitations in 
the ability to distinguish from each other the regular and 
chaotic patterns in a finite lattice. The identification 
of these patterns is further obscured by the fact that 
they may be superimposed on each other (as shown be- 
low). It should be stressed that the Peres' method is 
not quantitative — it does not yield (at least not directly) 
a calculable measure of quantum chaos which could be 
compared with other measures like, e.g., the Brody pa- 
rameter. 

In spite of these limitations, however, the method has 
a great potential to disclose important features of the 
mechanisms governing the breakdown of integrability and 
rise of chaos in low-dimensional systems. Its great ad- 
vantage is that the structural information on individ- 
ual eigenstates is represented by a single variable (the 
average of a suitably chosen observable), which allows 
one to use a simple visualization technique incorporat- 
ing simultaneously a large number of states. Note that 
in higher than 2D cases, the spectral lattice would have 
to be drawn in a multidimensional space, which would 
require to develop a sophisticated computer software for 
pattern recognition. Here, as we only deal with two- 
dimensional systems, the most efficient software is that 
already built in the human brain. 

Peres originally introduced his method in a more gen- 
eral way. He started from the simple fact that the time 
average of an arbitrary classical observable is a trivial in- 
tegral of motion (similarly as any function in the phase 
space which assigns a constant value to all points of the 
same trajectory). This makes it possible, for an arbi- 
trary system, to create an unlimited number of integrals 
of motion. Of course, this does not alter the fact that 
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FIG. 1: Peres lattices for J = eigenstates of the geometric collective model in the "2D even" quantization (see Sec. III). 
The points represent individual eigenstates with coordinates Ei (energy) and {P)i (expectation value of the respective Peres 
operator). Two Peres operators, and H' , are employed, the results shown in the left and right columns, respectively. Row 
(a) corresponds to the fully integrable case, 5 = 0. Rows (b)-(d) depict the disturbance of the lattice with gradually increasing 
nonintegrable perturbation until reaching the most chaotic case, B — 0.24 (e). All quantities are given in relative units. 



the system is nonintegrable, in general. Indeed, the func- 
tions corresponding to the new integrals are singular in 
the chaotic part of the phase space, hence they do not 
generally allow one to construct a transformation to the 
action-angle variables. Nevertheless, the dependence of 
time averages represents an interesting probe into the 
system's dynamics at given energy. 



The time averaging can be applied in quantum me- 
chanics, as well. Let us take an arbitrary Hermitian 
operator P, which in the present context will be called 
Peres operator. One can construct an operator P asso- 
ciated with the time average of quantity P. This oper- 
ator has the property that the time-averaged expecta- 
tion value (P)|^(o)) = limT^oo ^ Jo {i^{'f^)\P\i^{t))dt for 



4 



an arbitrary initial state |'^(0)) can be calculated as 
{P)mo)) = (V^(0)|P|V^(0)). The operator P is readily 
obtained from the Heisenberg image P}i{t) of P through 

P= \\m I- f P^{t)dt. (2) 

It is now straightforward to see that P fulfills the com- 
mutation relations [i^, P] =0, hence is an integral of 
motion. 

Peres showed that in an integrable system the set of 
points Ei versus Pi (where Pi is a fixed value of P in the 
zth eigenstate) forms a smoothly deformed regular lattice. 
This is so irrespective of the choice of the operator P 
used for evaluating the averages. Even if we choose an 
observable P which is not an integral of motion, P ^ 
the corresponding lattice for an integrable system will be 
ordered. The proof makes use of the EBK quantization 
and the fact that in an integrable system any additional 
integral of motion (including P) must be a function of the 
actions (for an integrable system, P is constant on the 
phase-space tori). Therefore, any distortion of regularity 
of the lattice signals the onset of chaotic motions. 

The expressions Pi and {P)i yield the same values and 
can be interchanged. We call these values P-averages, 
while the set of points Ei versus {P)i for an arbitrary (in- 
tegrable or nonintegrable) Hamiltonian is denoted here as 
the Peres lattice. Note that in Refs.[7, 9j, a more pictorial 
term "quantum web" was proposed. 

We want to stress that there is no restriction in the 
choice of the Peres operator P. Different choices give 
different lattices, but the separation of levels into reg- 
ular and chaotic parts of the lattice is independent of 
the choice. This consequence of the Percival conjecture 
will be discussed below. It is illustrated in Fig. [l] where 
the right-hand column shows lattices for another Peres 
operator than that used in the left-hand column. The 
rows correspond to the same values of the perturbation 
strength. We observe that the overall degree of chaos in 
each adjacent pair of images is about the same. More- 
over, it can be shown (cf. Fig.[2| that the states allocated 
in the regular (chaotic) part of one lattice lie in the reg- 
ular (chaotic) part of the other lattice, as well. [Note 
that illusive differences in the numbers of points on both 
sides of Fig. [l] (and some of the forthcoming figures) are 
caused by eventual accumulation of multiple points with 
very close coordinates.] 

III. HAMILTONIAN AND PERES OPERATORS 

In this section we briefly introduce the Hamiltonian of 
the Geometric Collective Model (GCM) of atomic nuclei 
[2^ in a nonrotating regime and suitable Peres operators. 
The model has been discussed in Part I [1 . The GCM 
Hamiltonian H = Hq + BH^ consists of the integrable 
part 

Ho = T-p^^p^ (3) 



and a nonintegrable perturbation 

i^' = /5^cos37, (4) 

where polar coordinates /3 and 7 stand for dynamical 
variables (shape parameters) and T is the kinetic en- 
ergy, which involves the associated momenta [1]. The 
corresponding Cartesian coordinates read as (x, y) = 
(/3cos7, /?sin7), with the momenta transformed accord- 
ingly. Parameter B is the perturbation strength, a 
model-specific version of the above- introduced general 
variable A. This choice corresponds to the GCM potential 
V = A(3^^B(3^cos3-f^C(3^ [1^ with (A,C) = (-!,+!). 
As C is positive, the Hamiltonian for any energy E de- 
scribes motions confined within a finite domain of p. The 
three degenerate global minima of the potential V are 
located at /3 > 0, 7 = |^ or (for P > or < 0, respec- 
tively), and a single local maximum is at = 0. As in 
Part I, all quantities are considered dimensionless. 

We take into account two different and physically rel- 
evant quantization schemes, which are connected with 
two- and five-dimensional versions of the system (here- 
after referred to as 2D and 5D cases, respectively) pi. 
The kinetic term T is different for both schemes, nev- 
ertheless in both cases it is proportional to the squared 
Planck constant over where K stands for an effective 
mass parameter of the system. The fraction = /K 
is called the classicality parameter. The value of this pa- 
rameter adjusts the absolute density of quantum spectra. 
In the following we set iiT = 1 and vary the value of h. 
The diagonalization of the Hamiltonian is performed in 
the appropriate 2D or 5D harmonic oscillator bases. The 
2D case is further split to even and odd case, referring 
to the symmetry or antisymmetry with respect to the 
7^—7 inversion. 

A set of eigenenergies and eigenvectors is obtained, for 
which the Peres lattice is constructed. We consider two 
types of the Peres operator. The first one is identified 
with the square of the angular momentum operator L 
connected with the rotations varying angle 7. In the 2D 
case, this is the Casimir invariant of the 0(2) algebra of 
rotations in the (/?, 7) plane: 

^2D = ^^^ n^m^, m = 0,3,6,..., (5) 

cf. Eq. (7) of Ref. [1 . Note that the eigenvalues, 
indicated in the last equation, involve multiples of 3 due 
to the required symmetry of the eigenfunctions with re- 
spect to rotations about angle 27r/3 [1^. In the 5D case, 
is the Casimir invariant of the GCM algebra 0(5) [15] 
restricted to value J = of the 0(3) angular momentum 
(null rotations in the ordinary space). We have 

^5D = ^^^sin37^ ^ n^v(v + 3), (6) 
sm 37 07 (77 

V = 0,3,6, 

cf. Eq. (5) of Ref. [1]. The eigenvalues of this operator 
are enumerated by an integer v (in the nuclear context 
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called seniority), which for J = again has only the 
values equal to multiples of 3 [15 . 

For the Hamiltonian eigenstates, and L^j^ take 
the values from Eqs. ([5| and ([6|, respectively, only for 
5 = 0, in which case the operators commute with the 
Hamiltonian. For 5 7^ 0, the energy eigenstates 7) 
are mixtures of states with different values of m or 'u. In 
any case, the average {L'^)i quantifies oscillations of the 
wave functions in the direction of angle 7 (examples 
shown below). 

The second Peres operator used in our analysis is iden- 
tified with the perturbation H' from Eq. Q. It is worth 
noting that the expectation value of H' in an eigenstate 
IV^i) of the general Hamiltonian ([T]) coincides with the 
derivative dEi/dX of the ith energy level at the given 
value of the control parameter. In the present case, the 
angular part of the Hq eigenstates has the property that 
(cos37)i = 0, hence 



_ dEi 



= 



(7) 



B=0 



(see the upper right panel of Fig.jT]). The vanishing slope 
at 5 = is consistent with the symmetry of the spectrum 
Ei{B) under the refiection B —B. If 5 7^ 0, however, 
the average is generally nonzero and satisfies {H')i = 
{Ei — {Ho)i)/B. The disturbances in the P = H' lattice 
therefore show up as departures of individual points from 
the line {H')i = 0, as seen in the right-hand column of 
Fig. [1] This facilitates the visual inspection of results. 



IV. RESULTS AND DISCUSSION 
A. Comparison of Peres lattices 

The effects accompanying the decay of regularity in 
Peres lattices were discussed from a general viewpoint 
in Sec. [Hj Figure [l} which corresponds to the GCM in 
the 2D even quantization with the above two choices of 
the Peres operator, illustrates the gradual transition from 
ordered to disordered lattices. Note that the case B = 
0.24 depicted in row (e) corresponds to the minimum of 
the classical regular fraction /^eg, see Part I [1 . The 
gradual onset of chaos in these lattices will be analyzed 
in more details in Sec. IIV CI 

In Fig. [2] we compare Peres lattices obtained (for both 
Peres operators) in the 2D even and 5D quantizations. 
The value B = 0.62 belongs to the island of strongly 
pronounced regularity close to the resonance of /3 and 7 



vibrations, see Sec. IV D The regularity shows up as a 
large area of ordered points, which starts at the lowest 
negative energies and spreads over to positive energies, 
where it is joined by a rising chaotic area. Despite the 
spectra for different quantizations show significant differ- 
ences [P, the form of Peres lattices is rather similar. 

In order to demonstrate the coincidence of the regular 
and chaotic regions in the lattices for different Peres op- 
erators, we have highlighted three of the states in the first 



row of Fig. [2j marking them by a square, a bullet, and 
a diamond. Probability densities for the corresponding 
wave functions are depicted in Fig. [3j The wave func- 
tions as well as the location of the respective points in 
the Peres lattice show that the square and the diamond 
correspond to regular states, while the bullet represents 
a chaotic state. In Fig. |2|a) we see that this assignment 
is consistent for both choices of the Peres operator. Let 
us note that both regular levels (a) and (c) in Fig. [s] ex- 
hibit a large increase of the wave-function magnitude in 
a region where a certain periodic trajectory oscillates in 
the classical case p^ . 

Although, as emphasized above, we can choose an ar- 
bitrary operator for plotting the Peres lattice. Fig. |2] in- 
dicates that some operators may be more suitable than 
others. For some choices, a part of the regular region 
in the lattice can pervade into the chaotic area and hide 
there behind a disordered mesh of points. In such cases, 
one cannot decide whether a level inside a chaotic re- 
gion is indeed chaotic. (On the other hand, overlapping 
regular areas form a regular area again.) While there is 
no doubt that these observations demonstrate limitations 
of the Peres method, one can improve its resolution by 
employing several incompatible Peres operators. Indeed, 
Fig. [2] (right) shows that for the three highlighted states 
a better choice of Peres operator is P = H' . 

Figure [4] displays Peres lattices for 2D even quantiza- 
tion with different classicality constant K:. In Part I, we 
have shown that by tuning the value of one scales the 
absolute density of quantum states but does not affect 
statistical properties of the spectra. Here we want to 
show that these changes do not infiuence the main fea- 
tures of the Peres lattice. The variations of the lattice 
with n = 25, 100 • 10~^ for both Peres operators are 
observed in rows (a)-(c) of Fig. [4j Note that B = 1.09 is 
a value of the control parameter for which the system ex- 
hibits very rich structures with well pronounced minima 
and maxima of the dependence classical regular fraction 
/reg on energy (see the insets). A decrease of the value 
increases the density of states (the system gets closer to 
the classical limit) and serves as a zoom into the sea of 
levels: one can see finer details of the lattice but (because 
of computational limits) a smaller fraction of the spec- 
trum is available. For a comparison, the box in all panels 
of the same column encloses a fixed region of energy x 
P-average. It is seen that the structures observed in the 
lattice become wealthier in details as k decreases, but 
the overall appearance of the relevant part of the lattice 
remains the same. 



B. Links to classical dynamics 

In Part I, we investigated the connection between the 
classical measure of regularity f^eg (the fraction of the 
regular phase-space volume) and the quantum measure 
represented by (1 — a;), where co stands for the Brody pa- 
rameter. Our conclusion was that both measures entail 
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FIG. 2: Peres lattices of the GCM at B = 0.62 with averages of (left) and H' (right) for 2D even (row a) and 5D (row b) 
quantizations {h = 5 - 10~^). In row (a), three states denoted by full symbols are identified in both lattices, demonstrating that 
the assignment of a given state to a regular or chaotic part of the lattice does not depend on the choice of the Peres operator. 
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FIG. 3: Squared wave functions for the three marked states from Fig. [2] (2D even quantization). The states (a) and (c), which 
correspond to the square and the diamond, respectively, are taken from the regular part of the lattice (the 1995th and 1885th 
level, respectively) and exhibit well pronounced quantum scar effects. The state (b), associated with the bullet, belongs to the 
chaotic part (the 1890th level) and shows an ergodic behavior. 



qualitatively the same energy-dependent behavior, irre- 
spective of the method of quantization. In the present 
context, new questions appear, namely: How strong is 
the correspondence between the behavior of f^eg and the 
character of Peres lattices? Or more specifically, is there 
a correlation between ordered (disordered) parts of the 
lattice and regular (chaotic) parts of the phase space? 

In order to find an answer we plot figures showing the 
dependence fj-eg{E) in the insets of Fig. |4] Pure visual 
inspection discloses strong correlations between the in- 
crease of /reg and the occurrence of regular domains in 
Peres lattices for both operators (L^) and {H'). 

Following the dependence in row (a) of Fig. [4] we find 
/reg = 1 for negative energies and observe fully regular 
patterns in the corresponding part of both lattices. This 



is the domain where the quadratic-well approximation 
of the GCM potential is valid. Note that the regular 
pattern at low energies is present in Peres lattices for all 
values of parameter B. At energies just below E = the 
classical regularity begins to drop (forming a small fold 
dit E ^ 0)^ which is manifested in the lattice of (L^) by 
a band of disordered points with a small tail penetrating 
to the regular area at zero energy. Passing through the 
totally chaotic area with /reg = 0atl<£^<3, the 
regularity begins to rise again. This is accompanied by a 
formation of a new regular pattern in the lattices at large 
values of (L^) or small values of {H'). 

One can switch to rows (b) and (c) of Fig. [I] and con- 
tinue in the same manner. A remarkable phenomenon 
appears at E ^ 35, where /reg reaches for a while the 
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FIG. 4: The GCM Peres lattices for B = 1.09 and different values of the Planck constant: h = 0.02 (a), 0.05 (b), and 0.1 
(c). The 2D even quantization was employed with the same Peres operators as in the previous figures. Lower values of k 
yield denser spectra, which can therefore be evaluated only in narrower energy intervals (numerical limitations). The box (red 
online) encloses the same area in the three panels of each column. Insets show the dependence of the classical regular fraction 
/reg on energy (see Part I T) for the energy domain displayed in each row. One can observe that the /reg dependences are 
correlated with ordered and disordered areas in the corresponding lattices in both columns. 



value of full regularity. This somewhat surprising behav- 
ior was discovered in Ref. [3 . Here we may trace the 
signatures of regularity in the two Peres lattices. For in- 
stance, the lattice in the lower right panel gets locally 
contracted to a narrow interval of {H') and develops a 
highly organized pattern for £^ > 35. Distortions of this 
pattern start appearing at energies E > 50, where the 
classical regularity is decreasing again. Let us note that 
for energies much above the range shown in Fig. |4j the 
order increases to the asymptotic value /reg 1 for all 
values of parameter B. This is due to the term of the 
potential which becomes increasingly important at high 
energies, generating predominantly regular dynamics. 

As explained in Sec. |n[ Peres invariants can be intro- 
duced on both classical and quantum levels. One can 
determine the classical analogue of the Peres operator P 
and calculate its average {P)c over an arbitrary trajec- 



tory. In this way, a function in the classical phase space 
can be constructed for any Peres invariant. In the right- 
hand panel of Fig. [sj the function {L'^)c is shown (coded 
in shades of gray) for a certain values of the control pa- 
rameter and energy on the y = section of the phase 
space. Remind that (x, y) represent Cartesian counter- 
parts of the polar coordinates (/?, 7), while (px^Py) are 
the corresponding momenta, and that there is no differ- 
ence between 2D and 5D cases on the classical level [l]. 

On the left-hand side of Fig. [5) the standard Poincare 
section is plotted for 100 crossing trajectories. We ob- 
serve that the chaotic area identified in the Poincare sec- 
tion is covered by one shade of gray in the map of {L'^)c- 
This follows from the ergodicity of chaotic motions, which 
ensures that any vicinity of each point in a chaotic phase- 
space domain is visited by a single trajectory. Therefore, 
almost the whole domain yields a single value of the clas- 
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FIG. 5: The y — section of the classical phase space (co- 
ordinate X versus momentum px) for B = 0.62 and E = 0.2. 
Left: Poincare section formed by crossings of 10^ randomly 
chosen trajectories with the plane of the section (10^ crossings 
for each trajectory). Right: Density plot of the Peres invari- 
ant (L^) calculated classically across the section (on a mesh 
of 500 X 500 points). Dark and light regions (low and high 
values of the Peres invariant, respectively) contain regular tra- 
jectories (see the left panel) and simultaneously correspond to 
regular domains in the Peres lattice [see Fig. [2] (a) dXE ^ 0.2, 
where the regular domains are located at the lower and upper 
sides of the lattice] . 



sical Peres invariant (with exceptions including periodic 
orbits that however fill only a zero-measure subset of the 
phase space). On the other hand, in the regular islands of 
the Poincare section the shade of the {L'^)c image grad- 
ually changes. 

The corresponding Peres lattice was shown in 
Fig. [2] (a), with squared wave functions of the selected 
states depicted in Fig. |3] The agreement with classical 
results in Fig. [5] is remarkable. The trajectory respon- 
sible for the "scar" of the wave function in Fig. [s] (a) 
passes the central regular part of the phase-space section 
in Fig. [5] close to {x^Px) = (—0-8,0), yielding a medium 
value of {L'^)c- Indeed, the respective state (denoted by 
the square) is localized in the medium part of the (L^) 
Peres lattice. The trajectory contributing to the wave 
function in Fig. [s] (c) falls to the dark regular regions of 
the density plot in Fig. [5) which again corresponds to the 
value of (L^) for the respective point (the diamond). A 
comparison of Figs, [s] and [2] (a) indicates an excellent cor- 
respondence between the results based on the Poincare 
and Peres methods (if the latter one is supplemented by a 
classical calculation as in the right-hand panel of Fig. [5| . 



C. Decay of regularity 

We return now to Fig. [l] in order to discuss in more de- 
tail the mechanism of the transition from the integrable 
dynamics at 5 = to the chaotic B > {) regime. As 
pointed out above, in the integrable case (row a) the 
quantity LF' is an integral of motion and {H') is identi- 
cally zero. 

Let us look at the curved chains of points apparent in 
panel (a) of Fig. [l] (left). These chains, which begin at 
(L^) = and lead upwards, connect states with a con- 



stant sum N = n(3-\- m^, where n/3 is the radial quantum 
number and m = 3m^ represents the angular-momentum 
quantum number (the latter increases towards the up- 
per end of the chain). If the perturbation is turned on, 
some of the points within the same chain start moving 
against each other, forming a kind of "condensation cen- 
ters"; see the left panels (b) and (c). A detailed inspec- 
tion discloses that the most affected levels lie in the short 
stretches of the chains which are nearly parallel with the 
vertical axis. In other words, the perturbation is most 
efficient for the states which are very close in energy and 
in the value of (i^^). Indeed, exploiting the perturbation 
theory, we can say that if the perturbation matrix ele- 
ment is nonzero (at 5 = 0, the operator couples only 
the states differing in by ±1), the proximity of levels 
leads to an increased mixing. This in the present case 
shows up as an attraction to a common value of (L^) for 
the whole bunch. For B = the vertical stretches are 
developed in the chains located within the energy inter- 
val ^ E < OA and, consequently, the corresponding 
levels in these chains are most vulnerable if B starts to 
increase. This is why small perturbations affect first only 
a very limited part of the lattice, as observed in panels 
(b) and (c). 

If we continue increasing the perturbation, more and 
more levels become influenced by the interaction. For 
sufficiently large values of the levels start interact- 
ing between neighboring chains and the whole structure 
gradually breaks down, see Fig. [1] (d). For B = 0.24 
(panel e) the lattice is totally disintegrated. We have 
just reached the most chaotic parameter region, where 
only the deepest levels form a regular lattice due to the 
validity of the quadratic-well approximation. 

The size of the perturbation can be quantified with the 
aid of the other Peres operator, i.e. by the value of {H')^ 
which is displayed in the right-hand column of Fig. [l] 
Rows (b) and (c) help to discover that not only the lev- 
els with ^ > 0, but also also a few of those with E < 
become disturbed by a small perturbation (this was not 
visible in the left-hand panels). For £^ > 0, we observe 
several regular arcs of points at {H') > and some more 
disordered points with {H') < 0. Both these groups of 
points correspond to the "condensation centers" appar- 
ent in the left-hand panel. The {H^) > part of the 
lattice contains states with (cos37)i > 0, hence 7 cen- 
tered around values 0, ^ and ^ (saddle points of the 
potential). On the other hand, the {H') < part col- 
lects states with (cos37)i < 0, hence 7 ~ |,7r, ^ (global 
minima of the potential). Examples of both these types 
of wave functions will be given later in Fig. [6] It is some- 
what surprising that the more regular part of the lattice 
is connected with the states localized in the saddle-point 
regions, whereas the states localized around the minima 
seem to be more chaotic. 

For moderate perturbation strengths, great majority of 
points in Fig. T] (right) remains located at (H^) = 0, in- 
dicating the absence of structural changes. These points 
correspond to the unperturbed parts of the lattice in the 
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FIG. 6: Squared wave functions of four successive eigenstates 
(states no. i = 52, . . . , 55 at E from 0.059 to 0.076) of the 
integrable B — Hamiltonian (for h — 0.02) in the 2D even 
quantization (the first row) and the same states for B — 0.005, 
0.01, and 0.05 (the second, third, and fourth rows, respec- 
tively). The respective values of both P-averages are given 
in each panel. The first three states (arranged in columns) 
are more sensitive to the perturbation than the fourth one, in 
accord with the lattices in Fig. [l] (see the text). 



left-hand column. As B increases, however, both posi- 
tive and negative halves of the {H') lattice become in- 
creasingly populated and finally the negative (irregular) 
part captivates absolute majority of points (row e). This 
agrees with the disordered form of the lattice in the left- 
hand column. 

The changes of the Hamiltonian eigenstates accompa- 
nying the above- described evolution of the Peres lattices 
are illustrated in Fig. |6] Its first row presents four un- 
perturbed {B = 0) wave functions (probability distri- 
butions in the 2D even quantization), while the second, 
third, and fourth rows demonstrate the effects of pertur- 
bation (for B = 0.005, 0.01 and 0.05, respectively) on the 
same states. The states correspond to four successive en- 
ergy levels, the associated values of both P-averages being 
given in each case. Note that the value of hi was chosen 
differently than above, so the states cannot be directly 
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FIG. 7: The lower and higher bounds (solid and dashed lines, 
respectively) of the classical Peres average {L^)c for ^ = 0. 
The bounds almost touch each other at B — 0.24, where 
the most chaotic region is located. Fluctuations of the lower 
bound near £^ = 0.5 are caused by the appearance and disap- 
pearance of small unstable tori. 



marked in Fig. [T] 

The rightmost column of Fig. |6] represents a state 
which is originally far away from the condensation cen- 
ters in Fig. [l] (left). Indeed, this state resists the smaller 
perturbation rather well. On the other hand, the most 
pronounced structural changes at 5 = 0.005 and 0.01 
are observed for the states in the first three columns of 
Fig. [6) The states in the middle two columns belong di- 
rectly to the condensation center of strongly interacting 
levels and the state in the leftmost column is close to 
it. One clearly observes the breakdown of the rotational 
symmetry and a gradual crossover to a trifoliolate form 
of the wave functions, particularly for the two states in 
the middle. While the states with {H') < are local- 
ized more around the minima of the potential, the ones 
with (H^) > dwell more in the saddle-point regions. In 
the fourth row of Fig. |6j which corresponds to the irreg- 
ular lattice at B = 0.05, cf. Fig. [l] (c), all four states 
are already perturbed. We observe that the form of the 
rightmost state has been transmitted to the third state 
in the last row and vice versa, as results from an avoided 
crossing of both levels. 

Finally, it is instructive to look also at the changes of 
classical Peres invariants with parameter B. We have 
calculated the lower and upper bounds of the classical 
average {L'^)c and show the results for £^ = in Fig. [?[ 
As we see, the interval of {L'^)c is contracted almost to 
a single value in the most chaotic case, B = 0.24, where 
the ergodicity is maximal. Surprisingly, even the small 
remnants of tori present there have the same value of 
(L^)c. If we step over the most chaotic point, several 
new tori with higher values of (L^) appear, which results 
in a widening of the interval between the bounds. Note 
that the narrow "neck" in (L^) at £^ = can also be 
observed in the quantum Peres lattice in Fig. [l](e). 
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D. A quasiregular region 

Looking at the form of (L^) and {H') lattices in row 
(e) of Fig. [l] one can notice two characteristic features: 
(i) The centroid value of (L^) exhibits roughly a linear 
increase with energy E. (ii) The lattice for {H') grows 
linearly only at low energies, while for higher energies 
it scatters around a roughly constant average. These 
types of dependences in both lattices are qualitatively 
understandable and remain approximately valid for all 
increasingly high values of parameter B. Nevertheless, 
the distribution of regular and irregular parts within the 
lattices of the above forms exhibits a high degree of vari- 
ability. 

At the first sight, one could expect that the increase of 
the perturbation strength in the GCM Hamiltonian from 
\B\ = {) io \B\ > {) should lead to a monotonous progres- 
sion of disorder, as described in the preceding subsection. 
Although this scenario is typical for many related systems 
of type ([T]), see e.g. Refs. [171 UHl^ it does not apply in 
the present case. 

The GCM is peculiar in two respects. First, as shown 
in Ref. [3] , the classical regular fraction /reg for negative 
energies converges to unity for asymptotically large val- 
ues of 5, when the GCM Hamiltonian can be rescaled to 
the form H^q = T + + cos 87. Let us stress that 
the type of order observed for £^ < in the B ^ 00 
limit is totally different from the B = case and that 
in the asymptotic limit the regularity fades away at pos- 
itive energies. Second, the competition between regular 
and chaotic dynamical modes gets surprisingly complex 
at medium values of \B\. The most important change in 
this range takes place around B ^ 0.6, where extensive 
regular patterns are established in the Peres lattices at 
both low and medium energies, see Fig. [2] (Note that 
the ordered dynamics at very high energies is connected 
with the dominance of the term of the potential [3 .) 

The B ^ 0.6 quasiregular region was briefly mentioned 
already in Part I [1 . As shown in Ref. [19 , it is closely 
related to a so-called "arc of regularity" observed [20] 
in the parameter space of the interacting boson model 
(IBM) [21 . Although the IBM is a more sophisticated 
model of nuclear collectivity that the GCM, it makes use 
of a similar language (involving quadrupole degrees of 
freedom) and yields comparable results. Since its discov- 
ery, the IBM arc of regularity has been subject to several 
analyses [19, [20l [22j [23]. Some interesting hints have 
been disclosed, but many questions remained open. 

In the GCM case, we observe a phenomenon very simi- 
lar to the IBM arc. A comparison of Figs.[l](e) and [2] (a) 
provides a clear evidence for a large increase of regular- 
ity between the two values of B. In fact, the pattern of 
ordered points, which dominates in the low-energy part 
of the lattice at B ^ 0.6, starts rising already before the 
maximum of regularity is reached and persists long after 
it is left. As an example, we show in Fig. ^the (L^) lat- 
tices for B = 0.52 (panel a) and B = 0.78 (panel b). It is 
obvious that the low-energy parts of both lattices exhibit 
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FIG. 8: Peres lattice for (L^) in 2D even quantization {h — 
5 • 10"^) at B = 0.52 (a) and 0.78 (b). We see that the 
regular pattern is present before and after the main peak of 
regularity dX B — 0.62, cf. Fig. |2](a). Both panels contain 
the same number of states. The insets depict selected wave 
functions (diamonds in the respective lattices) demonstrating 
the presence of [3 and 7 vibrations (only the sector around 
the minimum 7 = 47r/3 is shown in both cases). The state 
in the upper panel (i = 17) is at ^ = —0.449, the lower one 
{i = 1292) at ^ = 0.021. 



a great deal of similarity with Fig. [2] (a). 

The mechanism behind the E < {) pattern of ordered 
points visible in all (L^) lattices at medium and large 
values oi\B\ is connected with a competition of two types 
of vibrations. To show this, we apply the quadratic-well 
approximation, valid for |5| > at low energies above 
the potential minimum. It relies on the local use of a 2D 
oscillator potential 



k[3 



(8) 



where and 70 stand for a position of the potential 

minimum, and /c/3 = {^^^ and = (^^^^ for the 

rigidity of the oscillator in (5 and 7 directions. 

It turns out that the horizontal chains of points with 
increasing energy, which can be observed in the lattices 
in Figs. [2] (a) and [s] (a,b), correspond to states with a 
growing number of vibration quanta n/3. The vertical 
arrangement of these chains, on the other hand, follows 
an increasing number of 7-vibration quanta n^. Such an 
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interplay of vibrational modes in both f3 and 7 directions 
represents the basic organization principle for the low- 
energy part of the (L^) lattices for \B\ > 0. Remind 
that this is essentially different from the B = situation, 
when the lattice was determined by vibrational modes 
in P and rotational modes in 7. An example of a /3 x 7 
vibrational state (its squared wave function) is shown in 
the inset of Fig. [s] (a). 

A simple calculation shows that at 5 = 2/3 = 0.66 one 
gets kp = = 12. We encounter a resonance of the local 
oscillator frequencies in /3 and 7 directions, which leads to 
an additional regularization of the lattice. Interestingly, 
due to mutual interactions between levels the degeneracy 
in a wide interval of energies becomes maximal already at 
B = 0.62, where the main peak of regularity takes place. 
The resonance is responsible for the "condensation" of 
the (L^) lattice at £^ < along a nearly horizontal line 
of multiple points apparent in Fig. [2j Although in Fig. [l] 
a similar phenomenon was linked to initiating the first 
seeds of disorder, its role in the present case is rather 
opposite: it helps to clean up some disarranged parts of 
the lattice. Let us stress that the proximity of levels im- 
plies rapid structural changes with no immediate relation 
to chaos. It can indicate a crossover to chaos as well as 



emergence of order (imagine the scenario from Sec. IV C 
played in the reverse direction — with B decreasing to 0). 

It needs to be stressed that at 5 ~ 0.6 the patterns 
emerging in the Peres lattice and in wave functions go 
far beyond the quadratic-well approximation. Indeed, 
as shown in Fig. [2j the ordered part of the lattice ex- 
ceeds to medium energies, where the approximation de- 
teriorates and even becomes completely invalid (this is 
certainly so at energy £^sad of a saddle point of the po- 
tential, where the three regions around the minima merge 
together and form a single connected area). Surprisingly, 
even at £^ > > £^sad a large fraction of states still keeps 
the form with well distinguished P and 7 vibrations. This 
is exemplified by a selected wave function in the inset of 
Fig.|8](b), where the vibrational pattern remains confined 
around the potential minimum despite the fact that the 
energetically accessible domain unifies all three sectors. 
Note that peculiar /3 x j vibrational structures connect- 
ing all sectors can be found in even higher eigenstates, cf. 
Fig.|3](c). One may therefore assume that the observed 
regular island at B ^ 0.6 is due to a fortunate coinci- 
dence of resonating /3 and 7 modes in both E < and 
E ^0 domains. 

Qualitatively the same explanation is valid also in the 
IBM. There, the degeneracy of /3 and 7 vibrations was 
noticed empirically [23^ and later supported by theoret- 
ical arguments [19]. A detailed analysis of the f3 and 
7 modes in the IBM framework is in preparation [24 . 
The present work provides an independent verification of 
this mechanism in the simpler GCM case. It also clearly 
manifests the influence (probably specific for the present 
form of potential) of the low-energy ordering of states on 
the spectrum at higher energies, which is significant for 
the large extension of the regular region. 



V. CONCLUSIONS 

In this paper, we have continued and exceeded the 
work presented in Part I [l], whose main purpose was 
to test the Bohigas conjecture for different quantization 
schemes under the condition of a strong variability of 
chaotic measures with energy. We have revitalized an 
almost 35 years old method by Peres [6 and showed its 
great potential in the field of quantum chaos and even 
beyond. 

Peres lattices provide an excellent viewpoint to the 
landscape around the border between classical and quan- 
tum physics. This is so especially for systems with two 
degrees of freedom whose lattices can be drawn as two- 
dimensional diagrams, in analogy with planar Poincare 
sections of such systems. If applied within the domain of 
quantum chaos, Peres' method enables one to distinguish 
regular and chaotic behaviors on the level of individual 
states or subsets of states within the same energy inter- 
val. This is in contrast to traditional methods based on 
spectral statistics that assign the same degree of chaos 
to all levels within the same interval. 

Quite naturally, there are some limitations of the 
method. As seen, regular and irregular parts of the lat- 
tice can in some cases be superimposed on each other, 
which hinders their correct resolution. In particular, the 
distinction of chaotic states may be ambiguous since a 
superposition of two or more regular patterns may seem 
irregular. Nevertheless, we showed that this problem can 
in principle be bypassed by constructing more lattices 
with different Peres operators. Their optimal choice, 
which unavoidably depends on the concrete system un- 
der consideration, should be subject to further study. Al- 
though the Peres' method does not directly yield a cal- 
culable measure of quantum chaos, it represents an im- 
portant indicator providing new insights into the origin 
of chaotic behavior. 

However, our intention in this paper was to go even 
beyond the scope of quantum chaos, demonstrating that 
Peres lattices represent an extremely efficient and eco- 
nomic tool for studying significant features in large en- 
sembles of eigenstates across the spectrum. Relevant 
properties of the wave functions can be read off from 
the expectation values of suitably chosen Peres opera- 
tors. Instead of analyzing each individual eigenstate and 
its wave function, one may look at the associated Peres 
lattice where the desired information is contained in a 
synoptical way. As an example, we were able to closely 
follow the breakdown of integrability of the system and 
the rise of a new type of order. We believe that the 
results presented here may encourage similar studies in 
other systems. 

The present work completes our long-term effort to 
map chaotic properties of the geometric collective model 
of nuclear physics [U |3|. A great advantage of the ge- 
ometric model (and also of the related interacting bo- 
son model) is the apparent conceptual simplicity encod- 
ing strikingly rich complexity of dynamics. Let us note 
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that the above simphfied models capture the main phe- 
nomenological features of nuclear collectivity which are 
presently beside a fully microscopic description. The 
study of disordered collective dynamics within these 
models may be considered as an attack to the problem of 
chaos in many-body systems from the direction perpen- 
dicular to the mean-field approach. 

An interactive survey of our main results can be found 
at the website [25l. 
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